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Abstract 

We perform a detailed study of semileptonic form-factors for the doubly heavy baryons in 
the framework of three-point NRQCD sum rules. The analysis of spin symmetry relations as 
well as numerical results on various exclusive decay modes of doubly heavy baryons are given. 

o : 
o ; 

^ 1 Introduction 

c : 

i— j | Besides an experimental search for an explanation of the phenomenon of electroweak symmetry 
^ 'breaking and new physics beyond the Standard Model, high energies and luminosities of current 
^ 'and future particle accelerators provide a possibility to observe rare processes with heavy quarks. 

.An interesting topic here is a study of physics of doubly heavy baryons. An analysis of dynamics 
£>■ !of these baryons can play a fundamental role in an extraction of primary parameters of weak 
^ | quark interactions. Due to distinctions between the QCD effects inside the doubly and singly 
CN ' heavy hadrons, one may strictly constrain incalculable nonperturbative quantities, entering different 
^ schemes of calculations. 

The real possibility of such experimental measurements was recently confirmed by CDF Col- 
1 laboration by the first observation of B c meson ffl. As predicted theoretically 0, this long-lived 
I state of b and c quarks has the production cross sections, mass and decay rates, which are com- 
i ; patible with the characteristic values for the doubly heavy hadrons. Thus, the experimental search 
^■for the doubly heavy baryons can also be successful. Of course, such the search would be more 
^ ■ strongly motivated if it would be supported by modern theoretical studies and evaluations of basic 
>• ! characteristics for the doubly heavy baryons. 

Some steps forward this program were already done. First, the production cross sections of 
^ ] doubly heavy baryons in hadron collisions at high energies of colliders and in fixed target experiments 
were calculated in the framework of perturbative QCD for the hard processes and factorization 
of soft term related to the nonperturbative binding of heavy quarks 0. Second, the lifetimes 
and branching fractions of some inclusive decay modes were evaluated in the Operator Product 
Expansion combined with the effective theory of heavy quarks, which results in series over the 
inverse heavy quark masses and relative velocities of heavy quarks inside the doubly heavy diquark 
II, Third, the families of doubly heavy baryons, which contain a set of narrow excited levels in 
addition to the ground state, were described in the framework of potential models 0. The picture 
of spectra, obtained in this analysis, is very similar to that of heavy quarkonia. Fourth, the QCD 
and NRQCD sum rules J7| were explored for the two-point baryonic currents in order to calculate 
the masses and couplings of doubly heavy baryons [§, §, [TOfl . And fifth, there are papers, where 
exclusive semileptonic and some nonleptonic decay modes of doubly heavy baryons in the framework 
of potential models and within the Bethe-Salpeter approach were analyzed [[□], |T^| . 

In the present paper we estimate form-factors for the semileptonic decays of doubly heavy 
baryons together with semileptonic and some nonleptonic decay modes in the framework of three- 
point NRQCD sum rules for the case of spin 1/2 - spin 1/2 - baryon transitions only. The estimates of 
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contributions due to spin 1/2-spin 3/2-baryon transitions are also given the basis of QCD superflavor 
symmetry, emerging in the limit of very large heavy quark masses. In the limit of zero recoil we 
derive the spin symmetry relations on the form-factors, governing the semileptonic transitions of 
doubly heavy baryons. The use of these symmetry relations greatly simplifies further evaluation of 
form-factors for doubly heavy baryon transitions. A detailed analysis of baryon couplings, needed to 
model the phenomenological part of three-point sum rules is also provided. Our further exposition 
of the obtained results is organized as follows. In Section 2 we discuss our choice of interpolating 
currents between vacuum and corresponding doubly heavy baryon state and give numerical estimates 
of baryonic couplings in the framework of two-point NRQCD sum rules. Section 3 is devoted to 
the description of the method used to calculate the form-factors of interest. Here we present a 
derivation of spin symmetry relations between form-factors in the limit of maximal invariant mass 
of leptonic pair and give analytical expressions for corresponding double spectral densities. In 
section 4 we present numerical estimates of the form-factors studied together with predictions for 
the semileptonic and some nonleptonic decay modes. And finally, section 5 contains our conclusion. 



2 Two point sum rules 

In this section we describe steps, required for the evaluation of baryonic constants, as they will be 
needed later to model the phenomenological part of three-point sum rules. The question, which 
should be solved first is the choice of the corresponding interpolating currents for the baryons under 
consideration. So, in the next subsection we discuss their various choices and comment on the 
merits of the prescription, used in this paper. 



2.1 Baryonic currents 

As was mentioned in Introduction, in this work we consider only spin 1/2 - spin 1/2 transitions of 
doubly heavy baryons. Hence, the discussion of baryonic currents later in this subsection will be 
restricted to this case. For baryons, containing two heavy quarks, there are two types of interpolating 
currents: 

1) The prescription with the explicit spin structure of the heavy diquark from the very beginning 
is given by 

J~'« QI = lQ lT CT l5 Qi'}q k e tjk , 

Jsfc = [Q iT Cr^ m Q j ] ■ 7 m 7 5 g fc £r,-fc, (1) 

Here C is the charge conjugation matrix with the properties CjJ^C^ 1 = — 7 M and C^C~ X = 75, 
i,j, k are color indices and r is a matrix in the flavor space. 

2) The currents, which require a further symmetrization of heavy diquark wave function, have the 
form 

J %Q , = : (QlC l5 q p )Q' 7 : (2) 

The currents of the second type can be easily related to those of first type by the Fierz rearrangement 
of quark fields and the further symmetrization or antisymmetrization of diquark wave-function, 



depending on the diquark spin state. The first type of these currents was considered in || [TO 
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while the currents of second type were discussed in ||. The evaluation of form-factors, describing 
semileptonic decays of doubly heavy baryons, what is the goal of this paper, is much easy if we 
use the currents of the second type and make the symmetrization or antisymmetrization of diquark 
wave-function at the end of calculation on the level of form-factors. This procedure will become clear 
in the section with our numerical results for the exclusive decay modes of doubly heavy baryons. The 
reasons here are the manifest symmetry relations for form-factors, discussed by us later and rather 
simple calculations of various spectral densities for three-point correlation functions, which can also 
be done in full QCD framework without performing complicated angular integrations. However, as 
was shown by the authors of ||, for the second type of currents it is difficult, in general, to achieve a 
stability of sum rules predictions for the both extracted mass and coupling of doubly heavy baryons. 
To the same time, all these difficulties are absent for the currents of first type. So, the conclusion, 
which one may do here, is that, the appropriate choice of baryonic currents depends on the problem, 
which you would like to solve. 

Thus, in view of the lack of experimental information on the masses of doubly heavy baryons 
we will use for them the results of two-point NRQCD sum rules for the first type of currents || TU 



Taking these mass values as input, we calculate further the baryonic couplings of the second type. 

Next, let us discuss the couplings of strange heavy baryons, as they are appear in the semilep- 
tonic decays of some doubly heavy baryons. The currents, describing these hadrons, also classified 
according to the symmetry properties of light diquark wave-function. There are two spin 1/2 A- type 
(antisymmetric in the q and s - quarksQ) and two E-type (symmetric in the q and s - quarks) HQET 
currents, namely 

Jm = (q T C l5 s)Q v , J A2 = (q T C 1 ° l5 s)Q v , (3) 
(A - type) 

= (gW^hSQ,, Jx2 = (q T C 1 ° 1 k s) 1 k l5 Q v , (4) 
(E - type) 

where Q v is the HQET heavy quark spinor, moving with velocity v. The baryons, described by these 
currents, belong to the same SU(4) multiplet, that have as its lowest level the J p = | + SU(3) octet. 
As, these hadrons contain only one heavy quark, they belong to the second level of the mentioned 
SU(4) multiplet. This level splits apart into two SU(3) multiplets, a 3, states of which Eq are 
antisymmetric under interchange of two light quarks and thus described by A - type currents, and 
6, states of which Eq are symmetric under interchange of light quark and described by E - type 
currents. Actually, there may be some mixing between the pure 3 and 6 states to form the physical 
Eq and Sq states^. So, in what follows we will not distinguish between Eq and Eg - baryons and 
will exploit the fact that both states have non- vanishing overlap with e al31 : {q^ l C , y^Qp)s 1 : current. 
In other words, what we suppose to calculate is the semileptonic branching ratio of some of doubly 
heavy baryons into both Eq and Eq - baryons. 

Now, let us briefly describe the two-point NRQCD sum rules, used for their evaluation. 



Here q denotes one of the light quarks u or d. 
2 They both have the same I, J and P quantum numbers. 
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2.2 Description of the method 

We start from the correlator of two baryonic currents with the half spin 

U(p 2 ) = i J d 4 xe^(0\T{J(x), J(0)}|0) = pT^p 2 ) + F 2 (p 2 ), (5) 

Performing the OPE expansion of this correlation function, we get a series, different terms of which 
give us the contributions of operators with various dimensions. So, for Fi (i, 2) functions we have 
the following expressions 

F*(p 2 ) = Fr rt (p 2 ) + ifVXw) + Ff (p 2 )(^G 2 ) + iTVK^g) + ... (6) 



To obtain theoretical expressions for the Wilson coefficients, standing in front of different operators, 
one typically uses the dispersion relation 

7T Jo W — t 

where t = k ■ v, p^ = + (mi + m 2 )v^ and p\ denotes the imaginary part in the physical region 
of corresponding Wilson coefficients in NRQCDf]. The calculation of spectral densities p\ proceeds 
through the use of Cutkosky rules |TJ| and for the case of e afBl : (Q^C^q^Qy : current and different 
quark masses was done in the QCD framework in 0. Here we use the results of this work. The 
needed NRQCD spectral densities were obtained by simple NRQCD expansion of corresponding 
QCD expressions and the results of this expansion could be found in Appendix A. 

To relate the NRQCD correlators to hadrons, we use the dispersion representation for the two- 
point function with the physical spectral density, given by appropriate resonance and continuum 
part. The coupling constants of doubly heavy baryons are defined by the following expression 

(0\J H \H(p)) = iZ H u{v, M H )e ip - x , (8) 

where p = M H v and the spinor field with four-velocity v and mass M H satisfies the equation 
f/>u(v, Mjj) = u(v,M H ). 

We suppose that the continuum part, starting from the threshold value w cont , is equal to that 
of calculated in the framework of NRQCD. Then, equalizing the correlators, calculated in NRQCD 
and given by the physical states, the integrations above w cont cancel each other in both sides of 
sum rules relation. Further, we write down the correlators at the deep under-threshold point 
to = —(mi + m 2 ) + t at t — > 0. 

Introducing the following notation for the n-th moment of two-point correlation function 

71" Jo [w + m 1 + m 2 ) n 
and using the approximation of single bound state pole, we can write the following relation 

Mi = l^J^-r (10) 



3 Here mi and are the heavy quark masses and v, L is four-velocity of the baryon under consideration. 
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From which one can read off the corresponding expression for the baryon coupling in the moment 
scheme 



I 7"W 1 2 
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where we see the dependence of sum rules on the scheme parameter n. Therefore, we will tend to 
find the region of parameter values, where, first, the result is stable under the variation of moment 
number n, and, second, the both correlators F\ and F 2 reproduce equal values of coupling constants. 
In the QCD sum rule analysis of || there was a problem, that values of baryon coupling constants, 
obtained from F\ and F 2 correlation functions significantly differ. To cure this problem, in || it 
was proposed to include in calculations also contributions, coming from the OPE expansion for the 
correlator of two quark fields 
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With an account of these corrections the quark condensate contribution to moments gets modified 
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The derivation of two-point HQET sum rules for the heavy baryons with the strangeness follows the 
same lines as that for baryons with two heavy quarks. Here we only comment on the differences. 
To obtain the HQET expressions for the spectral densities we again use the QCD result of ||. 
However, the transition between the QCD expressions for the doubly heavy baryons and HQET 
expressions for the heavy baryons with the strangeness is more intricate. First, we should take a 
limit when one of the heavy quark masses goes to zero and second, we should allow this quark to 
condense. We have calculated explicitly the s -quark condensate contribution to the both F\ and 
F 2 correlation functions and subtracted l/m s poles from gluon condensate contribution, related to 
the strange quark condensate due to the following heavy quark expansion 
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The appearing logarithmic singularities can be related to the mixed quark condensate with the help 
of the same heavy quark expansion 
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And, finally, one must subtract the nonsingular gluon condensate contribution, belonging to the 
quark condensate, what can be easily done with the calculated explicit expression for quark con- 
densate contribution. The resulted HQET spectral densities for the case of ordinary baryons with 
strangeness were collected by us in Appendix A. For the strange quark condensate contribution, as 
in the case of light quark condensate, we also take into account the corrections, coming from the 
OPE expansion for the correlator of two strange quark fields 
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With this corrections the s-quark contribution to the moments for the F\ and F 2 correlation func- 
tions has the form 

Mf(n) = -]m^-^^M s %n + l) + ^(^-m 2 \m s { -^^M Ss (n + 3), (17) 



Ml s {n) = M Ss (n) 



n\ 48 V 4 

ml - 2m 2 s (n + 2) 



16 



M Ss (n + 2) + 



(vr2 (^)_| m 2 K _ m 2 ))(n + 4)! 



where 



and 
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(w + m Q + m s ) n ' 
(ss)(m s + w) 2 {2rriQ + m s + w) 2 



(18) 
(19) 
(20) 



47r(m Q + m s + w) 

Now, having all theoretical expressions for baryon coupling constants in the moment scheme, we 
will proceed in the next subsection with the numerical estimates. 
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Figure 1: a) The values of \Zn\ 2 and \Z V n\ 2 H£ 6 - baryon couplings as functions of the moment 

number n; b) the value of 1 - 1 2 ' " 1 average as function of the moment number n for the case of 
2£ b - baryons. 
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2.3 Numerical estimates 

In this subsection we present the results on the coupling constants of doubly heavy baryons and the 
heavy baryons with the strangeness. In the scheme of moments, which we employ here to extract 
the baryonic couplings from the two-point sum rules, the dominant uncertainty in estimates comes 
from the variation of heavy quark mass values. In the analysis we chose the following region of 
quark mass values 

m b = 4.6 - 4.7 GeV, m c = 1.35 - 1.40 GeV, (21) 
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Figure 2: The results for the S£ c - baryon interpolating current in the case of decaying 6-quark (the 
(c T C7 5 g)6-current ): a) the values of \Zn\ 2 and \Zn\ 2 H£ c - baryon couplings as functions of the 

moment number n; b) the value of " 1 2 " 1 average as function of the moment number n for the 
case of El c - baryons. 
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Figure 3: The results for the S£ c - baryon interpolating current in the case of decaying c-quark (the 
(6 T C7 5 g)c-current ): a) the values of \Zn\ 2 and \Zn\ 2 H£ c - baryon couplings as functions of the 

moment number n; b) the value of ]{Zn • ■ average as function of the moment number n for the 
case of El c - baryons. 



what is the ordinary choice used in sum rules estimates of heavy quarkonia. For the strange quark 
mass we use the value m s = 0.15 GeV. 

Next point, we would like to discuss, is an account of Coulomb corrections inside the doubly 
heavy diquark. As is well known, these corrections give large contribution to baryon coupling 
constants and are essential for relative contributions of perturbative and condensate terms to the 



correlator With a good accuracy at low or moderate values of moment number, the effect of 
Coulomb interactions can be written as overall Zommerfeld factor in front of perturbative spectral 
density of heavy subsystem for the square of baryon coupling. But, as it was shown in the same 
Zommerfeld factors should be taken into account in calculations of three-point correlation functions 
considered later in this paper. It occurs, that for the form-factors for the semileptonic transitions of 
doubly heavy baryons these corrections cancel each other in average. So, the calculation of desired 
form-factors for the doubly heavy baryons can be consistently performed without accounting for 
Coulomb corrections either, provided we neglect them both in the two-point and three-point sum 
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Figure 4: The values of \Zn\ 2 and \Zn\ 2 S* c - baryon couplings as functions of the moment number 
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Figure 5: a) The values of \Zn\ 2 and \Zn\ 2 S£ - baryon couplings as functions of the moment 

number n; b) the value of ^ Zn ■ ^ Zn ■ average as function of the moment number n for the case of 
S£ - baryons. 



rules. This is the approach, we will follow in the present work for the evaluation of form-factors. 

The dependence of estimates on the threshold of continuum contribution in the two-point sum 
rules is not so valuable as on quark masses. We fix the region of w cont as 



w cont = 1.3-1.4 GeV, 



(22) 



which is in agreement with our previous estimates of doubly heavy baryon coupling of the first type 
currents in the same framework of two-point NRQCD sum rules. For the condensates of quark and 
gluons we use the following regions: 



(qq) = -(250 - 270 MeV) 3 , m 2 = 0.75 - 0.85 GeV , (—G ) = (1.5 - 2) • 10" 2 GeV 4 , (23) 

7T 



and 



(ss) = 0.8 ± 0.2 {qq) 



(24) 
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Figure 6: a) The values of |zj^| 2 and \Zn\ 2 - baryon couplings as functions of the moment 

mber n; b) 
- baryons. 



number n; b) the value of ^ Zn • ^ Zn ■ average as function of the moment number n for the case of 



As we already mentioned, for the second type of currents used here, we evaluate the coupling 
constants only and use the masses of doubly heavy baryons, calculated by us previously 0, as 
inputs 



Mu = 3.47 ± 0.05 GeV, M= h = 6.80 ± 0.05 GeV, 



10.07 ±0.09 GeV, 



(25) 



which are in agreement with the values obtained in the framework of potential models. For the 
masses of heavy baryons with the strangeness, appearing as products of semileptonic decays of some 
of the doubly heavy baryons, we use the following values: 



5.8 GeV, M Sc = 2.45 GeV. 



(26) 



Figs. 1-6 show the dependence of baryon couplings on the momentum number. We find that the 
stability regions for l-Z^)] 2 determined from the F\ and F 2 correlators coincide with those, obtained 
in the analysis of two-point correlation functions for the first type of currents || . However, for some 
of the couplings, calculated here, we see a sizeable difference in the predictions coming from the F\ 
and F2 correlation functions. This problem could not be alleviated by the variation of parameters. 
So, in order to determine the corresponding coupling values to be used in the phenomenological 
part of three-point sum rules, we consider an average coupling for these currents, whose square is 
given by the average of squares for Z 1 and Z 2 couplings. The resulted values of baryonic coupling 
are 
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Having estimates for the couplings of initial and final state baryons with respect to semileptonic 
transitions, we will continue in the next section with the determination of form-factors. 
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3 Three-point sum rules 



In this section we describe our framework for the calculation of form-factors, governing the semilep- 
tonic decays of doubly heavy baryons. Here we derive the spin symmetry relations between various 
form-factors, arising in the limit of the maximal invariant mass of leptonic pair and give analytical 
expressions for corresponding spectral densities. 




Figure 7: The diagram, corresponding to the three-point correlation function considered in the 
paper. 

Following the standard procedure for the evaluation of form-factors in the framework of QCD 
sum rules, we consider the three-point function 

U„ = i 2 j d 4 xd 4 y(0\T{J HF (x)UO)J Hl }\0)e^ x e- vp 'y, (29) 

where J M is the vector or axial transition current, matrix elements of which between baryonic ground 
states we would like to calculate. Fig. 7 shows a diagram, corresponding to the mentioned three- 
point function. The theoretical expression for the three-point correlation function can be easily 
calculated with the use of double dispertion relation 

nf eor \ Sl , s 2 , q 2 ) = — !— f°° dsi [°° ds 2 . Bt^^ll} + subtractions, (30) 

where the desired spectral density could be obtained with the help of Cutkosky rules [TJ]]. We 
will continue with the calculation of spectral densities later in this paper after discussing the spin 
symmetry relations for form-factors. The latter, as will be seen, greatly simplify the calculations to 
be done. Now, let us discuss the phenomenological part of three point sum rules under consideration. 
Saturating the channels of initial and final state hadrons by ground states of corresponding baryons, 
we have the following phenomenological expression for the three-point correlation function: 



(Hj( Pi )\J Hi \0) 



s 2 — Mir 

spins ^ n F 



(H F (p F )\J u \Hj( Pl )}- 



4 - Mf Jl 
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The formfactors for the weak spin ~ - spin | baryon transitions are usually modeled as following: 

(HpMlJ^Hzipj)) = (32) 
u(p F ) [ llt (F? + Ff-ft) + ia^F? + F 2 A 75 ) + q^Fl + F 3 A l5 )]u( Pl ) 

However, in the NRQCD limit it is more convenient to use an alternative parametrization 



(HpMlJ^Hjipj)) = (33) 
u(p F ) ( 7M Gr + vfcl + < Gl + 75(7^ + vffii + «; G3 )) u(p,), 



fl Z 1 /i d 1 |U\ //* 1 fit Z 1 

where these two parametrizations can be related to each other with the help of the following relations: 

F?(t) = G\ + {m F + mi ) (-^G v 2 +-^—Gl 



2m j 2rriF 



Ff(t) = -Gf - (m F - mj) G# + G% 



— - -i- 

2m/ 2m/ 

L cr + ^ 

2m/ 2m/r 

2m/ 2mp 

^(*) = ——G2 + —^—Gf , 
2m/ 2m p 

Naively, all these six formfactors in either parametrization are independent, but, as we will show in 
the next subsection, in the limit of zero recoil the semileptonic decays of doubly heavy baryons can 
be described by the only universal function, an analogue of Isgur-Wise function. 



3.1 Symmetry relations 

Now, let us discuss the spin symmetry relations among the form-factors, arising in the limit of zero 
recoil for the final state baryon. That is, we consider a limitQ, where vj 7^ Vp and u = (vj ■ v F ) — > 1. 
The theoretical expression for the three-point correlation function for the case of heavy to heavy 
underlying quark transition in this limit has the following form: 



where 



n gfc«r) „ C IW (q 2 )(l + fp) (l _ ^ + ^ (35) 



vj = vi + — -(vj - v F ) (36) 
2mi 

v F = V p + ^-{v F - vi). (37) 



4 For the discussion of this limit see [151 
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So, for this type of transitions we have already, from the very beginning, the only universal func- 
tion and no further analysis is required. The theoretical expression for the three-point correlation 
function in the case of heavy to light underlying quark transition has more complicated form 

nf-) ~ {^(g 2 )^ + Uq 2 )^f + Uq 2 )h,a - 7 5 )(1 + h) (38) 

Considering different convolutions of the theoretical and phenomenological three-point correlation 
functions with Lorenz structures made of hadron velocities and 7 - matrices and equating them, we 
obtain two relations on the semileptonic form-factors in this case 

[GX + GX + GX) = e W {q 2 ) (39) 
Gf = £™V) (40) 

and a relation between £i(g 2 ) functions 

Recalling also, that in any considered transition we always have heavy baryons in initial and final 
states and requiring that, the appropriate projections do not change the theoretical expression for 
the three-point correlation function, we may conclude, that in this limit there are only two form- 
factors of order unity G\ = Gf = £(w), while all others are suppressed by heavy quark masses. 

Having derived the spin symmetry relations, we came to situation where we should calculate 
the only universal function in order to obtain estimates on semileptonic or nonleptonic transitions 
of doubly heavy baryons. 



3.2 Spectral densities 

As, we said before, the calculation of spectral densities is straightforward with the use of Cutkosky 
rules for the quark propagators. However, the resulted expressions for NRQCD spectral densities 
are different^ for the cases of heavy to heavy or heavy to light underlying quark transitions, so below 
we have classified the calculated quantities. For the trace of correlation function with v we have 

1) heavy to heavy underlying transition 



pP ert = 1 2 =(m\-Am\ml + ?>ml - Amly/sJ + Smxmly/s} + 

VA(s/, s F ,q 2 ) 

6m 2 S/ - 4m 2 s 7 - 4m lS 3 / 2 + sj + Am\ log ^~ m% \ (41) 

m 3 

p „ = 4_p^ 

(mi + m 3 )^\(s I ,s F ,q 2 ) 

2) heavy to light underlying transition 

5 It is simply an artifact of NRQCD approximation. 
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4(2tt) 2 v / S jA(sj, s f , ? 2 ) 

g 2 + 2m 2 \/^7 + s f) ~ Bm^m! — v /sJ) 2 (2m 2 + m 2 — 2g 2 + 4m 2 y / i7 + (43) 
2s F ) + ml(9ml + 2m\ - 9q 2 + 18m 2v /s7 + 9s F ) + \2m\(m\ - q 2 + 
/JJ — mi 



2m 2 ^/s] + s F ) log 



m 3 



P qq = ~- ^^===(ml-q 2 + 2m 2y ^+s F -ml)(qq), (44) 

(mi +m 3 )y/X(s I ,s F ,q 2 ) 



where 



X(x, y, z) = x 2 + y 2 + z 2 - 2xy - 2xz - 2yz. (45) 
and the integration region in the double dispertion relation is determined by the condition 

-1 < 1 = ((s f + s F - g 2 )(s 7 + m 2 - m\) - 2 Sl (s F - m\ + m 2 )) < I. 

V A(sj, sf, g 2 )A(s/, mf , m§) 

The notations in the above expressions should be clear from Fig. 7. Having derived theoretical 
expressions for the three-point correlation function, we may proceed now with the evaluation of 
form-factors. In numerical estimates we will use the Borel scheme for the form-factor extraction and 
so, below we give an expression determining the universal Isgur-Wise function for the semileptonic 
decays of doubly heavy baryons 

1 i rs? r^ F h 

^ IW ^ 2 ) = (n^2o M M 7 7 / p{si,S F ,q 2 )d Sl ds F X 



'(mi+mj) 2 J {m\+m,2) 2 

si-M 2 s F -M 2 F , 
- w —) exp( w 



exp(-^ 2 -^) exp(-^-^), (46) 



where Bj and B F are the Borel parameters in the initial and final state channels. 



4 Numerical results 

In this section we give the results of numerical estimates on the form- factors for the spin 1/2 - 
spin 1/2 doubly heavy baryon transitions. Assuming the pole resonance model for the dependence 
of mentioned form-factors on the square of lepton pair momentum we make predictions on the 
semileptonic, pion and p - meson decay modes. 

4.1 Form-factors 

The analysis of NRQCD sum rules in the Borel scheme gives us the estimates of the value of Isgur- 
Wise (IW) function at zero recoil for different types of transitions between doubly heavy baryons, 
shown in Table 1. 

In Figs. 9-11 we have plotted the dependence of the normalization of IW-function on the Borel 
parameters in the channels of initial and final state baryons. Exploring the stability of NRQCD 



13 



sum rules upon variation of these parameters just give us the results quoted above. The subtle 
point in the presented analysis is the choice of the threshold values in the baryon channels. In the 
present analysis we put the same values as in the analysis of two-point sum rules. The results on 
the formfactors and later their comparison with the results of potential models convince us, that 
we made a right choice. However, the latter are in general different from the ones used in two-point 
sum rule analysis in moment scheme. For the situation, where it is the case we refer the reader to 
PI- 




Figure 8: The value of £(1) for the transition E^ b — > S£ c as function of Borel parameters in the sj 
and sp channels. 

For the sake of comparison, we also provide here the estimates of the values of IW-function 
at zero recoil performed by us in the framework of potential models (PM). In this approach the 
normalization of IW-function is given by the overlap of initial and final state baryon wave-functions. 
For simplicity we assume the factorization of doubly heavy baryon wave-function in the diquark and 
light quark - diquark wave-functions. From the HQET-description of heavy mesons we know, that 
the light quark affects the normalization of IW-function at zero recoil only in 1/toq order. Thus, 
in our case its affect can be neglected and the normalization of IW-function is given by the overlap 
of diquark wave-functions. Taking the Gaussian anzaz for the diquark wave-function we get 



2W X Wy 

vol + vol 



V 



«"(!)= I ^S) • < 47 > 



where 



and 



Z/| 2 V /3 f\Z F \ 2 ^ 1/3 



Wx = 2lT \i2 ' w y = 27T ( ~it ) (4 ^ 



tPM\ 



2^1^(0)^(0)1- (49) 

Here we have related the parameters of diquark wavefunctions to the baryon couplings, obtained 
previously by us in the framework of two-point NRQCD sum rules. In table 1 we have gathered the 
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Figure 9: The value of £(1) for the transition S£ c — > as function of Borel parameters in the sj 
and sp channels. 

results of sum rules and PM on the values of IW-functions at zero recoil. We see, that within the 
errors of sum rules method (15 %) the obtained results are very close to those of PM. 



Mode 


£(1) SR 


£(1) PM 




0.85 


0.91 


1 — 'be ' 1 — >cc 


0.91 


0.99 


1 — be ' 1 — >bs 


0.9 


0.99 


' — CC ' 1 — >cs 


0.99 


1. 



Table 1: The normalization of Isgur-Wise function for different baryon transitions at zero recoil. 

Next, to obtain the dependence of formfactors on the square of momentum transfer we exploit 
the pole resonance model. So, for the IW-function we have the following expression: 

£™V)=£o-^, (50) 
1 — — i — 

m poU 

with 

n^poie = 6.3 GeV for the b — > c transitions 
wipc-ie — 1-85 GeV for the c — > s transitions. 

4.2 Semileptonic decays 

Now knowing all formfactors, describing semileptonic transitions of doubly heavy baryons we can 
estimate the semileptonic decay ratios for the transitions under consideration 

/W m ax ^"P 
d^^QQ'^^QQ'), (51) 
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Figure 10: The value of £(1) for the transition H£ c H£ as function of Borel parameters in the sj 
and sf channels. 



where 



«^ = ^nr^; q 2 = M I+ M F -2M I M F w. (52) 



2. \ /,.!/,. 



For the ^P- we have 

<toj 



where 



o?w du> ' 



- ^W^W + WI, (54) 

S<^- ^ Sfc) = ^W^^^d^l" + (55, 

Here H\ Fi \ w = H XpXw — x , where and \w are helicities of final state baryon and W - 
boson correspondingly and the functions obey the following symmetry relations: 

The functions remained after the application of this relation can be further expressed in terms of 
calculated in previous subsection IW - functions with the help of the following formulae 

Hi/£,i = -2^/MrMA^V)i IW {w) (57) 

H$ >0 = -^=^2MjM F (w T l)(Mj ± M F )£ IW {w) (58) 



The results of numerical estimates, done with the help of presented formulae can be found in Table 
2. 
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Figure 11: The value of £(1) for the transition as function of Borel parameters in the sj 

and sf channels. 



To calculate the ir or p - meson decays we assume the hypothesis of factorization [M. The 
corresponding formulae for the decays of doubly heavy baryons with a pion or p - meson in the 
final state can be easily obtained from those for the semileptonic decays by a simple substitution of 
leptonic tensor by the it - meson current tensor f%pZp* or p - meson current tensor f^m? p {—g^v + 



m 



2\ 



(( Mj + M F f - q 2 ) dT 
(Mj + Mp) 2 dq 2 " 



F Hl ^H F n = 67T f^ifi) — -J-2\q 2 =rnl, (59) 



^^H FP = -^0^{(M I -M F ) 2 ((M I + M F ) 2 -q 2 ) + 

dT 



2m%(M I - Mp) 2 - q 2 )}^\ q2=m> (60) 

where cii(/f) = ^j-(C + (p)(N c + 1) + C-(ji)(N c — 1)) and N c = 3 is the number of colors. In 
numerical calculations we put a\ = 1.26. The results for these nonleptonic transitions can be also 
found in Table 2. To calculate the branching ratios for exclusive decay modes we used the values of 
doubly heavy baryon lifetimes, calculated by us previously |J. There is some difference in concrete 
numerical values of lifetimes, obtained in different papers j|, ||. In || we have commented on the 
uncertainties in the resulting values of lifetimes related to the heavy quark mass values. There 
is, however, one more uncertainty remained, connected with the value of light quark - diquark 
wave-function at origin. In present there are two approaches to estimate this value: 1) assuming, 
that this value is the same as the value of D-meson wave-function at origin; 2) extracting this 
value from the comparison of hyper-fine splittings in doubly and singly heavy baryons. Here we 
used the estimates for the lifetimes made in the second approach, as they are the most complete 
ones. The values, presented in Table 2 already include the contribution of spin 1/2-spin 3/2 decay 
channels. To estimate the latter we have used the results of |TT[, where the contribution of these 
channels was calculated for the case of — > S cc + W baryon transition, and assumed, that, 
according to superflavor symmetry, it constitutes 30 % from the contribution of corresponding spin 
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1/2-spin 1/2 transitions for all transitions between doubly heavy baryons. In calculations of H£ 6 
and - baryon decay modes we have taken into account a factor 2 due to Pauli principle for 
the identical heavy quarks in the initial channel. In the case of S£ c — > S^ c X-baryon transition the 
same factor comes from the positive Pauli interference of the c-quark, being a product of 6-quark 
decay, with the c-quark from the initial baryon. Here, we also would like to mention, that for the 
Hb c -baryon decays the mentioned positive Pauli interference contribution is dominant among the 
other nonspectator contributions^, so we do not introduce other corrections here. However, in the 
case of S+ + — > X- baryon transition the negative Pauli interference plays the dominant role and 
thus should be accounted for explicitly. From the previously done OPE analysis for doubly heavy 
baryon lifetimes f|, [| we conclude that the corresponding correction factor in this case is 0.62. We 
would like also give a small comment on our notations. The Hq s in Table 2 stays for the sum of Hq 
and Sg decay channels. 



Mode 


Br (%) 


Mode 


Br (%) 


—bb ~^ —bc lu l 


14.9 


^te "> 


4.9 


__. w+ir, 

—be —cc lu l 


4.6 


-be ^ -bs iV l 


4.4 


-be E bs lu l 


4.1 


s+ + - s+> 


16.8 


Zte - ^cM 


7.5 


—bb ~* —bc n 


2.2 


—bb ~^ —bcP 


5.7 


' — be * ' — 'cc ^ 


0.7 


1 — 'be * 1 — 'ec^ 


0.7 


—be * —cc P 


1.9 


' — 'be * ' — ccP 


1.7 


—be —bs n 


7.7 


—be ~^ —bs n 


7.1 


-be —bsP 


21.7 


—be ^ —bsP 


20.1 


' — 'cc * ' — 'cs^ 


15.7 


5+ _^ H° 7T+ 

cc cs 


11.2 


' — 'cc * ' — csP 


46.8 


' — 'cc * ' — csP 


33.6 







Table 2: Branching ratios for the different decay modes of doubly heavy baryons. 



The previous studies of exclusive decays of doubly heavy baryons JTT], l2]| exploited the spin-flavor 
symmetry of QCD ]T7|], arising at the very large quark mass limit. The fundamental representation 
of the SU(6)®U(1) spin-flavor symmetry group consists of two-component heavy quark spinor, 
scalar and axial-vector di-antiquark fields. This representation can be explicitly written in terms of 
nine-component vector with the four- velocity v. 




where A% satisfies the constraint v = and the effective lagrangian for this field is 

Ceff = g ^vM iv- < ~D > * w 



'Here we use the results of OPE analysis for the inclusive decay modes of doubly heavy baryons [Q, || 
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where Ai is a 9x9 mass matrixQ, given by the following expression: 



M 





and 



_ 7 ° 

%, = I 10 



OOg 

Here (7 = diag(l,-l,-l,-l) is the usual metric tensor. Next, to make connection with the hadronic 
states, one considers a tensor product of \fr v with one light antiquark field. Thus, this hadronic 
super mult iplet puts together singly heavy mesons and doubly heavy antibaryons. However, such 
supermultiplet is not completely flavor-independent even in the heavy quark mass limit, as there 
remains internal mass-dependent heavy di-antiquark dynamics. The singly heavy baryons with the 
strangeness in the discussed approach belong to the different supermultiplet and this fact should 
be taken into account, when calculating form-factors for the semileptonic transitions between dou- 
bly heavy and singly heavy baryons. Such analysis within the framework of potential models was 
performed previously by M.A.Sanchis-Lozano [IT] for the case of 5£ c -baryon decays. There, to cal- 



culate the form-factors, it was assumed that the latter are given by the overlap of Coulomb diquark 
wave-functions with small non-perturbative corrections, given by the presence of a light quark in 
the baryons under consideration. It is just the approach we have used in our PM estimates. To 
reduce the number of independent form-factors, there was performed an analysis of spin-symmetry 
relations between various form-factors in the limit of zero recoil. The author, using different argu- 
ments, had came to the same conclusion as we have did in the present work. That is, all semileptonic 
transitions of doubly heavy baryons are governed by the only universal function, an analogue of 
Isgur-Wise function. The numerical results on the normalization of Isgur-Wise function at zero 
recoil completely agree with our estimates both in the framework of potential models and NRQCD 
sum rules. The given predictions for the semileptonic decay modes of 5£ b -baryons nicely agree with 
the ones presented in this paper, taking into account correction factor due to the wrong values 
of S£ c -baryon lifetime used in that paper. There is also a paper, where the diquark semileptonic 
transitions where calculated within the Bethe-Salpeter approach. However, the numerical results 
presented in this paper are very strange. It is suffice to say, that, for example, according to these 
results the semileptonic branching ratios of E^-baryon decays should be approximately 50 %, what 
is very unlikely. 

To finish the discussion of the obtained results we would like to note, that the latter are also in 
agreement with the estimates of inclusive decay channels performed by us previously H. H . 



5 Conclusion 

In this paper we have presented the analysis of exclusive decays of doubly heavy baryons in the 
framework of NRQCD sum rules. We have provided complete numerical study of baryonic couplings 
and semileptonic form-factors. The values of semileptonic and some nonleptonic exclusive modes 

7 Note, that the particular form of the mass matrix depends on the fields normalizations. 
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are also given. To conclude, we would like to discuss what also can be and should be done in 
the study of exclusive decays of doubly heavy baryons. First, it will be instructive to perform the 
similar analysis for the baryon currents of the first type. We have check, that these two schemes of 
calculation give the similar results only in the case of S£ b -baryon decays. Second, one may perform 
an analysis of doubly heavy baryon exclusive decays in full QCD and not relay on the pole resonance 
model for the form-factors. We plan to present the results of such analysis in nearest future. And, 
third one may try to calculate the lifetimes of doubly heavy baryons in the framework of QCD sum 
rules. It is a very interesting task, as we will explicitly see the effect of large nonspectator effects, 
studied previously in the OPE framework, on various exclusive modes. 
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Appendix A 

In this Appendix we have collected theoretical expressions for spectral densities of Wilson coeffi- 
cients, standing in front of various operators, obtained as the result of OPE expansion of two-point 
correlation function. 

For the case of e a/3A : (Q'JC^q^Qx : current we have 



pert 
Pi 



oert 



HI 

n 



p? 



pf 



pT 



mix 
Pi 



mix 
P2 



2\^2^mim 2 (mi + m 2 ) 
105(mi + m 2 ) 3 7r 3 



w 7 / 2 (m 1 m 2 (12m 2 — 13w) + hm\w + m 2 (12m 2 + 5w)) 



2\[2Jm\m,2(m\ + m 2 ) 7/2/ x 2 2 ,^ n „ NN 

— — — ; w 1 mim 2 12m 2 - w) + m 2 w + m 1 (12m 2 + bwj) 

105(mi + m 2 ) tt 6 . 



y / m 1 m 2 (m 1 + m 2 ) 

AV2(m 1 + m 2 ) 3 7r 

y / m l m 2 (m 1 + m 2 ) 

4\/2(m 1 + m 2 ) 2 7r 
1 



V / w(mim 2 (4m 2 — 5u>) + hm\w + m 2 (4m 2 + 

V / w(mim 2 (4m 2 — tu) + m 2 -u; + m 2 (4m 2 + it;)) 

V / w(2m2tf 2 — 28m^m 2 w(4m 2 + iu) + 



1536-N/27r(mim 2 ) 3/2 (mi + m 2 ) 7 / 2 
mim 2 it;(— 16m 2 + 35w) + ml(32m 2 l + 8m 2 w — u> 2 ) — 
8m 2 m 3 (8m 2 - 13m 2 w + 28w 2 ) + m 3 (-96m^ + 2YJm\w 2 )) 

' / w(122m2W 2 - 100m^m 2 (4m 2 + iu) + 



7680 v / 27r(mim 2 ) 3/2 (mi + m 2 ) 5 / 2 
mim 2 u>(880m 2 + 327u>) — 40m 1 m 2 (8m 2 — 41m 2 w — 20u> ) + 
5m 2 (32m 2 + 8m 2 tu — u> 2 ) + 5m 3 m 2 (— 96m 2 , + 64m 2 u> + 169w 2 )) 

(m\m\ (64m 2 - 397w) + 



(61) 
(62) 

(63) 

(64) 

(65) 



(66) 



2048v / 27r(mim 2 ) 1 /2( mi + m 2 ) 7 / 2 v^; 
30mim 2 (4m 2 - 21w) + 105m 2 w + m?(-104m 2 + 17iu) + 10m?m 2 (-16m 2 + 19iu)) 



(67) 



: (2m 3 m 2 (96m 2 — 43u>) — 



2048 v / 27r(m 1 m 2 ) 1 /2( mi + m 2 )5/2^ 
6m 1 m 2 s (4m 2 — 23w) + Ihm^w + m^(104m 2 — 17w) + m 2 m 2 (64m 2 + 45w)) 



(68) 



Here m\ is the mass of Q-quark and m 2 is the mass of Q'-quark. For the case of e al3X : (Q^C^q^Sx : 
current we have 



oert 



W 



•ert 



Pf 



135m s w — 12mQm s w(5m s + 3w) + 4mg(5m s + hm s w + w )) 
(9m s w 2 + 5mg(4m s + w) — rriQw(5m s + w)) 



80m 2 Q 7i 3 

m s w 3 
40mQ7r 3 

— - — 3 — (— 85m 2 -u7 3 — m 2 Qw{3m s + 2w) 2 + mQm s w 2 (33m s + 26w) + 
Q 

m^(ms + Am s w + 2w 2 )) 



(69) 
(70) 

(71) 
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m Q m s 



r (2mg(m s + w) + w 2 (m s + w) + mQw(2m s + 3w)) 



47r(mg + w) 2 

— T (-2185m 2 s w 3 - 28m 2 Q w(3m s + 2w) 2 + 2m Q m s w 2 (A37m s + 3Uw) - 
32rriQ(m s + w) + 32m 3 Q (m 2 s + Am s w + 2w 2 )) 
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8mg(m s + 4-u; — m s log 8) + 2Am 2 Qm s {rriQ + m s + w) log it; ) 
^-(mj + 168m> 2 - 5m 3 Q (m s + w) - 3m Q m s w{l8m s + 19w) + 

m 2 Q (10m 2 s + 22m s w + llw 2 )) 
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